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蛇 足 石 杉 COBRA 基因 家 族 的 分 子 生物 信息 学 及 表达 分 析 


HEW, EJ, RAE, ZEHN 

(广西 中 医药 大 学 药学 院 ， 南 宁 530200) 
JH 要 : FARER COBRA 基因 家 族 成 员 分 子 生物 信息 学 特征 及 组 织 表 达 规 律 ， 为 
COBRA 基因 进一步 研究 提供 依据 。 基 于 蛇 足 石 杉 的 全 长 转录 组 数据 ， 通 过 生物 信息 学 技术 
对 该 家 族 成 员 (HsCOBRAs) 的 理化 性 质 、 结 构 域 、 保 守 基 序 、 顺 式 作用 元 件 、 基 因 表 达 量 
等 进行 分 析 。 结 果 表 明 : (1) ERE Ae RAP SE i aE HH 24 个 HSCOBRAs 家 族 成 
员 ， 其 中 酸性 蛋白 9 个， 稳定 蛋白 I+, HKERAS 个， 具有 跨 膜 结构 的 蛋白 7 个 ， 具 
有 信号 肽 的 蛋白 3 个 。 (2) 亚 细 胞 定位 在 细胞 壁 、 叶 绿 体 、 细 胞 核 、 细 胞 膜 上 。 (3) 结构 
分 析 发 现 HsCOBRAs 有 7 种 结构 域 和 6 种 保守 基 序 ,部 分 成 员 具 有 高 度 保守 的 CCVS 结构 。 
(4) HsCOBRAs 具有 CAAT-box、TATA-box 等 46 种 顺 式 作用 元 件 。 (5) HsCOBRA2 在 叶 、 
孢子 、 茎 、 芽 胞 中 的 表达 量 均 最 高 。 蛇 足 石 杉 COBR4 基因 家 族 分 子 生 物 信 息 学 及 表达 特征 
可 为 HsCOBRAs 的 进一步 研究 及 生物 学 功能 验证 等 提供 理论 依据 。 
KRR: EEB, COBRA 基因 家 族 ， 生 物 信息 学 ， 全 长 转录 组 ， 表 达 分 析 
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COBRA gene family in Huperzia serrata 


HUANG Yumei, TENG Jianbei, TU Dongping', LIANG Liuguan 
(College of Pharmacy, Guangxi University of Traditional Chinese Medicine, Nanning 530200, 
China) 
Abstract: The paper clarifies the molecular bioinformatics characteristics and tissue expression 
patterns of the COBRA gene family members of Huperzia serrata, to provide a basis for further 
research on COBRA genes. Based on the full-length transcriptome data of the Huperzia serrata, 
the physicochemical properties, domains, conserved motifs, cis-acting elements, and genes 
expression of the family members ( HsCOBRAs ) were analyzed by bioinformatics techniques. 
The results were as follows: (1) A total of 24 HsCOBRAs family members were screened in the 
full-length transcriptome of Huperzia serrata, including 9 acidic proteins, 11 stabilizing proteins, 
5 hydrophobic proteins, 7 proteins with transmembrane structures, and 3 proteins with signal 
peptides. (2) Subcellular localization was found in the cell wall, chloroplast, nucleus, and cell 
membrane. (3) Structural analysis revealed that HsCOBRAs have 7 domains and 6 conserved 
motifs, and partial members have a highly conserved CCVS structure. (4) HsCOBRAs have 46 
cis-acting elements such as CAAT-box and TATA-box. (5) HsCOBRA2 had the highest 


expressionin in leaves, spores, stems and gemma. The molecular bioinformatics and expression 
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characterization of the COBRA gene family of Huperzia serrata can provide theoretical basis for 
further research and biological function verification of Hs COBRAs. 
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COBRA 基因 最 早 发 现 于 拟 南 芥 根 细胞 中 (Borner et al., 2003)， 是 一 类 广泛 存在 于 高 等 植 
物 和 藻类 植物 ( 任 昂 彦 等 ,2021) 中 的 纤维 素 合成 酶 基因 (Kasirajan et al., 2018)， 其 编码 的 蛋白 
盟 于 糖 基础 脂 酰 肌 醇 锚 定 蛋白 ( glycosyl-phosphatidyl inositol, GPD ， 包 含 N 端 信号 肽 、 碳 水 
化 合 物 结合 域 、 富 含 半 胱 氨 酸 的 CCVS 结构 域 和 C 末端 GPI 蛋白 的 @- 位 点 (Roudier et al., 
2002)。 目 前 ,学 者 们 已 鉴定 了 棉花 、 大 白菜 、 茶 树 、 番 前 、 玉 米 等 多 种 植物 的 COBRA 基因 。 
研究 表明 COBRA 可 通过 调节 纤维 素 微 纤 丝 的 合成 及 沉积 ,在 初生 和 次 生 细 胞 壁 纤维 素 的 生 
物 合成 中 发 挥 核心 作用 (Li etal., 2019; EMES, 2021), FF ZSIR XE. nb. HERR 
实 的 发 育 ， 影 响 着 植物 的 生物 量 (元 志 成 等 ，2020; Zaheer et al., 2022; Li et al., 2022)。 如 棉 
花 GhCOBL9A 的 过 量 表 达 可 使 生物 量 增加 (Niu et al., 2018); 茶树 CsCOBRA 基因 可 通过 调控 
叶片 细胞 壁 及 机 械 强度 ， 从 而 影响 茶叶 的 产量 ( 艾 安 涛 等 ，2021)。 
蛇 足 石 杉 (Huperzia serrata) PAMPER, ATR. ibd EL. i 
WES. TAPAS. BARETA APAT Bett. a AEE. BJ. 5e 
AE RERE, 2023) . HEBER A A Bal RY Bll A HELI Calzheimer’s disease, AD) 
RE BLUE IT AG HIAS RRM, 2007). IFS) 2050 年 ， 全 球 AD 患 病人 
数 将 达 1.39 亿 (Jia et al., 2018)， 我 国 老年 人 患 病人 数 将 达 3 003 (ERAS, 2019), WEE 
石 杉 需 求 量 不 断 增 大 ,市场 前 景 良好 。 近 年 来 ， 蛇 足 石 杉 的 人 工 栽 培 取 得 一 定 进 展 ,但 因 划 
生长 缓慢 ,繁殖 困难 ， 至 今 尚 未 实现 大 规模 栽培 ， 产量 远 远 不 能 满足 市 场 需求 。 高 产量 、 高 
效益 的 蛇 足 石 杉 现代 化 生产 已 成 迫切 需要 ,通过 COBRA 等 纤维 素 合成 酶 基因 的 调控 有 望 成 
为 蛇 足 石 杉 生 物产 量 增长 的 新 途径 。 
本 研究 以 蛇 足 石 杉 为 研究 材料 ， 从 其 全 长 转录 组 测序 数据 中 筛选 COBRA 基因 ， 通 过 生 
物 信 息 学 技术 对 基因 特征 及 在 组 织 中 的 表达 量 进 行 分 析 ， 探 讨 : (1) te ALAA COBRA 基 
因 家 族 的 生物 信息 学 特征 ; (2) COBRA 基因 在 蛇 足 石 杉 生长 发 育 中 可 能 发 挥 的 作用 ， 为 将 
来 蛇 足 石 杉 COBRA 基因 功能 的 鉴定 及 调控 葛 定 基础 ,为 蛇 足 石 杉 的 分 子 育种 从 而 提高 产量 
提供 依据 。 
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1 材料 与 方法 


1.1 材料 

蛇 足 石 杉 植株 于 2019 年 8 月 在 广西 南宁 市 大 明山 采集 ， 经 广西 药 用 植物 园 冯 世 讲 高 级 
工程 师 鉴 定 为 石 杉 科 植 物 蛇 足 石 杉 (Huperzia serrata) 的 全 草 。 蛇 足 石 杉 全 长 转录 组 数据 由 深 
圳 华 大 基因 股份 有 限 公 司 PacBio Sequel 平台 提供 。 
1.2 方法 
1.2.1 REA COBRA 基因 家 族 成 员 鉴定 

基于 蛇 足 石 杉 的 全 长 转录 组 测序 数据 及 其 注释 结果 ， 利 用 拟 南 芥 COBRA 基因 家 族 成 员 
的 DNA 序列 作为 种 子 序列 , BLAST 比 对 搜索 蛇 足 石 杉 全 基因 序列 数据 库 , de REL E 1075, 
3k fa RE LE ER COBRA 基因 家 族 成 员 。 
1.2.2 WERI COBRA 基因 家 族 的 理化 性 质 、 信 号 肽 、 亚 细胞 定位 和 结构 域 分 析 

将 蛇 足 石 杉 COBRA 基因 家 族 的 原始 碱 基 序列 通过 开放 阅读 框 ORF Finder 

(http://www.ncbi.nlm.nih.gov/orffinder) 获 得 潜在 的 蛋白 质 编 码 片段 。 利 用 ExPAsy(Expasy - 

ProtParam tool ) 对 COBRA 基因 编码 的 蛋白质 序 列 进 行 理 化 性 质 预 测 ，SignalP-6.0 


Iw 


C https://services.healthtech.dtu.dk/service.php?SignalP ) 进行 信号 肽 预测 ， Plant-mPLoc 
Chttp:/Avww.csbio.sjtu.edu.cn/bioinf/plant-multi/) 进行 亚 细 胞 定位 预测 ，Batch CD-Search 
C https:/Avww.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi ) 进行 结构 域 预测 ， 并 利 
TBtools1.112 软件 进行 可 视 化 处 理 。 
1.2.3 WEAF COBRA 基因 家 族 的 二 级 结构 、 三 级 结构 和 顺 式 作用 元 件 分 析 
分 别 通 过 Prabi(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)、 
phyre? (http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index ) 进行 二 级 结构 和 三 级 结构 
预测 。 以 蛇 足 石 杉 COBRA 家 族 基因 上 游 2 000 pb 序列 作为 启动 子 区 域 ， 利 用 Plant CARE 
Chttps://bioinformatics.psb.ugent.be/webtools/plantcare/html/) 进行 顺 式 作 用 元 件 预测 ， 算 选 
具有 注释 功能 的 顺 式 作用 元 件 ， 利 用 TBtools1.112 软件 进行 可 视 化 处 理 。 
1.2.4 te AZ COBRA 基因 家 族 的 保守 基 序 、 表 达 量 和 系统 发 育 分 析 
通过 MEME Suite 5.5.1 (https://meme-suite.org/meme/tools/meme ) 预测 基因 家 族 的 保守 
基 序 , 利用 TBtools1.112 软件 进行 可 视 化 处 理 。 通 过 计算 每 千 个 碱 基 的 转录 每 百 万 映射 读 取 
的 片段 (fragments per kilobase of exon model per million mapped fragments, FPKM) 均值 ， 分 
Nt COBRA 3&DS AO EE EAI. J. AL FFRAE, FAY TBtools 1.112 软件 绘 
制 基因 表达 量 热 图 并 进行 美化 。 通 过 MEGAT1 软件 对 有 蛋白 序列 进行 多 重 比较 ， 使 用 邻接 法 
(neighbor-joining) 构 建 系统 进化 树 (bootstrap=1 000)， 利 用 iTOL Chttp://itol.embl.de) 对 进 
化 树 进 行 美化 。 


2 结果 与 分 析 


2.1 HEA COBRA 基因 家 族 的 理化 性 质 分 析 

从 蛇 足 石 杉 全 长 转录 组 数据 中 共 第 选 出 24 条 COBRA 基因 家 族 的 碱 基 序 列 ， 分 别 命 名 
为 HsCOBRA1~HsCOBRA424， 其 碱 基 序列 可 通过 图 1: A 获得 。 基 于 蛇 足 石 杉 转录 组 测序 及 
ORFfinder 获得 氨基 酸 序列 (图 1: BO 。 利 用 ExPAsy 在 线 软 件 对 HsCOBRAs 编码 的 蛋白 质 
进行 理化 性 质 预测 ， 结 果 见 表 1。24 条 HsCOBRAs 蛋白 序列 中 ， 和 蛋白 氨基 酸 数量 为 73 
(HsCOBRA19) ~663 (HsCOBRAI4) 个 。 分 子 量 在 8 494.91 CHsCOBRAIS) ~73 084.23 
(HsCOBRA14) Da 之 间 ， 平 均 分 子 量 为 26 012.36。9 个 蛋白 (HsSCOBRA1、HsCOBRA5、 
HsCOBRA6 、HsCOBRA7 、HsCOBRA9 、HsCOBRA14 、HsCOBRA16 、HsCOBRA17 、 
HsCOBRA18 ) 等 电 点 在 4.53~6.97 之 间 , JJRETESR EH: ER 15 个 蛋白 的 等 电 点 在 8.17~10.13 
之 间 ， 为 碱 性 蛋白 。11 个 蛋白 (HsCOBRAI,. HsCOBRA2. HsCOBRA3., HsCOBRA4, 
HsCOBRA6 、HsCOBRA7 、HsCOBRA11 、HsCOBRA12 、HsCOBRA13 、HsCOBRA14、 
HsCOBRA23) 为 稳定 蛋白 ， 其 余 13 个 蛋白 均 为 不 稳定 蛋白 。HsCOBRAs 家 族 成 员 的 脂 溶 
性 指数 介 于 61.70~92.37 之 间 ; HsSCOBRA7、HsCOBRA8、HsCOBRA9、HsCOBRA12、 
HsCOBRAIS Jy ifi ZK TES A; 其 余 和 蛋白 为 杀 水 性 和 蛋白。 所 有 和 蛋白 均 带 有 负电 残 基 3 


(HsCOBRA19) ~77 (HsCOBRA16) 个 及 正 电 残 基 6 (HsCOBRA15) ~53 (HsCOBRAS ) 
个 
| 


o 


图 1 蛇 足 石 杉 COBRA 家 族 的 碱 基 (A) 及 氨基 酸 序列 B) 二 维 码 
Fig.1 The two-dimensional code of COBRA family base sequence (A) and amino acid 


sequence (B) of Huperzia serrata 


#1 HEH COBRA 家 族 理化 性 质 分 析 


Table 1 Physicochemical characteristics of COBRA family of Huperzia serrata 


氨基 酸 数量 分 子 量 PEER RERA RKE WAE WEB 
基因 名 Number of Molecular ”等 电 点 数 Aliphatic ”平均 系 残 基 残 基 
Gene name amino acids weight pl Instability index 数 (Asp + (Arg + 
(Da) index GRAVY Glu) Lys) 
HsCOBRAI 302 32 530.59 5.47 36.85 65.26 -0.420 37 28 
HsCOBRA2 471 52 393.40 8.17 38.96 81.80 -0.042 28 32 
HsCOBRA3 503 55 930.98 8.68 38.03 68.09 -0.246 32 42 
HsCOBRAA 471 52 393.40 8.17 38.96 81.80 -0.042 28 32 
HsCOBRAS 497 54 926.58 5.20 46.82 83.38 -0.208 66 53 
HsCOBRA6 122 13 210.99 4.99 32.01 74.43 -0.009 15 10 
HsCOBRAT 120 13 043.86 5.02 31.76 78.92 0.041 16 10 
HsCOBRAS 142 15 656.18 8.78 53.72 90.00 0.008 7 11 
HsCOBRA9 137 15 292.54 4.58 42.90 91.75 0.120 12 8 
HsCOBRA10 81 9 655.12 9.03 42.43 69.88 -0.267 7 10 
HsCOBRAII 255 28 653.78 824 39.64 68.47 -0.209 15 18 
HsCOBRA12 154 16 979.79 8.35 38.91 79.81 0.102 8 11 
HsCOBRA13 342 37 845.55 8.66 39.42 75.61 -0.080 17 24 
HsCOBRAI4 663 73 084.23 5.05 35.02 83.91 -0.076 57 43 
HsCOBRA15 75 8 494.91 8.74 56.17 71.47 0.128 4 6 
HsCOBRA16 80 9 176.60 6.52 53.67 80.50 -0.094 TI 7 
HsCOBRA17 127 13 913.37 4.53 47.61 71.34 -0.412 22 8 
HsCOBRA18 174 19 703.91 6.97 46.39 84.08 -0.010 12 12 
HsCOBRA19 73 8 531.23 10.13 61.70 80.27 -0.255 3 13 
HsCOBRA20 141 15 644.73 9.57 65.95 61.70 -0.397 7 16 
HsCOBRA21 119 13 578.72 9.52 42.76 65.55 -0.606 13 21 
HsCOBRA22 232 25 716.16 8.47 51.62 73.10 -0.656 29 31 
HsCOBRA23 194 22 823.72 9.23 18.98 92.37 -0.343 24 30 
HsCOBRA24 130 15 116.39 8.95 61.45 85.38 -0.189 6 11 


2.2 WEE AAS COBRA 基因 家 族 亚 细 胞 定位 及 信和 号 肽 分 析 
跨 膜 结构 预测 发 现 7 个 蛋白 (HsSCOBRA2、HsCOBRA3、HsCOBRA4、HsCOBRA9、 
HsCOBRA11、HsCOBRA12、HsCOBRA14) 分 别 具 有 1 个 跨 膜 结构 ， 其 余 均 无 跨 膜 结构 。 
言 号 肽 预测 发 现 HsSCOBRA2、HsCOBRA4、HsCOBRA14 具有 典型 的 GPI 锚 定 蛋 白 N 端 信 
号 肽 ， 信 号 肽 指数 分 别 为 0.975 8、0.975 8、0.999 8， 其 切割 位 点 分 别 可 能 在 N 端的 第 41 
和 第 42、 第 AT 和 第 42、 第 28 和 第 29 个 氨基 酸 残 基 间 ; 其 余 信 号 肽 指数 均 小 于 0.5. NEZH 
胞 定位 发 现 HsCOBRAs 家 族 成 员 分 布 广泛 ，1 个 (HsCOBRA9) 定位 于 细胞 壁 ，4 个 
(CHsCOBRA6、HsCOBRA7、HsCOBRA10、HsCOBRA23 ) 定 位 于 叶绿体 , 10 个 (HsCOBRA1、 
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表 2 MELA COBRA 家 族 亚 细 胞 定位 及 信和 号 肽 分 析 
Table 2 Subcellular localization and signal peptide analysis of COBRA family of Huperzia 
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蛋白 结构 域 分 析 结 果 显 示 HsCOBRAs XO 
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2) ， 分 别 为 COBRA、 


COBRA super family, PI31 Prot N, PI31 Prot C. PI31 Prot N super family. RING Ubox 


super family. TilS super family. 


其 中 COBRA 与 COBRA super family 分 别 分 布 在 4 个 蛋白 


序列 中 , 均 注 释 为 糖 磷脂 酰 肌 醇 。 除 HsCOBRAL 含有 2 个 结构 域外 , 13 个 蛋白 (HsCOBRA2、 
HsCOBRA3、HsCOBRA4、HsCOBRAS、HsCOBRA6、HsCOBRA7、HsCOBRA8、HsCOBRA12、 
HsCOBRA13、HsCOBRA14、HsCOBRA16、HsCOBRA17、HsCOBRA18) 仅 有 一 个 结构 域 ， 


其 余 的 10 个 蛋白 无 结构 域 。 
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图 2 蛇 足 石 杉 COBRA 家 族 结构 域 预测 图 
Fig. 2 The COBRA family domain prediction of Huperzia serrata 


2.4 he EA COBRA 基因 家 族 二 级 结构 分 析 
利用 Prabi 软件 对 蛇 足 石 杉 COBRA 基因 家 族 进 行 二 级 结构 预测 ， 表 3 结果 表明 
HsCOBRAs 蛋白 的 二 级 结构 均 由 无 规则 卷曲 、a 螺 旋 、 延 长 链 、 及 B 折 又 组 成 ， 所 占 比 平均 
值 分 别 为 46.02%. 29.3496. 18.95%, 5.70%. 6 个 蛋白 (HsCOBRA10、HsCOBRA19、 
HsCOBRA20 、HsCOBRA21 HsCOBRA22. HsCOBRA23) 中 a 螺旋 占 比 最 大 ( 占 
38.27%~58.25%) ， 其 次 为 无 规则 卷曲 ;其 余 HsCOBRAs 蛋白 中 的 无 规则 卷曲 比例 最 大 ( 占 
38.67%~64.29%) ; ATA RAM BA eA SIBI EG lee Crh 1.55%~14.81%) 。 
#3 We RURT COBRA 家 族 二 级 结构 分 析 
Table 3 The secondary structure analysis of COBRA family of Huperzia serrata 


pur 


蛋白 a 螺旋 延长 链 pris 无 规则 卷 
Protein Alpha helix (96) ^ Extended strand (96) Beta turn (96) | Random coil (%) 
HsCOBRAI 20.86 14.90 3.64 60.60 
HsCOBRA2 14.44 23.35 4.67 57.54 
HsCOBRA3 19.68 21.47 4.37 54.47 
HsCOBRA4 14.44 23.35 4.67 57.54 
HsCOBRA5 43.46 9.26 3.02 44.27 
HsCOBRA6 31.15 18.03 1.64 49.18 
HsCOBRA7 26.67 22.50 5.83 45.00 
HsCOBRA8 32.39 14.79 7.75 45.07 
HsCOBRA9 27.74 22.63 6.57 43.07 
HsCOBRA10 38.27 20.99 14.81 25.93 
HsCOBRAII 17.65 20.39 3.92 58.04 
HsCOBRA12 19.48 13.64 2.60 64.29 
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13.45 
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16.00 
26.25 
33.86 
25.29 
57.53 
43.97 
46.22 
48.71 
58.25 
12.31 


23.39 
22.62 
33.33 
22.50 
8.66 
21.26 
5.48 
17.02 
16.81 
10.34 
16.49 
31.54 


2.5 EAX COBRA 基因 家 族 三 级 结构 分 析 
利用 phyre? 软件 对 蛇 足 石 杉 COBRA 家 族 成 员 进行 三 级 结构 预测 ， 表 4 和 图 3 结果 表 


明 24 ME 


HsCOBRA16、HsCOBRA19、 


HsCOBRA7 均 具 


均 具 


HsCOBRA6、HsCOBRA7 4 


AKARE 


HsCOBRAS 序列 同 源 性 为 100%. 
表 4 EER COBRA 家 族 三 级 结构 分 析 
Table 4 The tertiary structure analysis of COBRA family of Huperzia serrata 


建 模 氨 基 酸 及 比例 
蛋白 Modeling amino 
Protein acids and 
proportions (%) 
HsCOBRAI 135(45) 
HsCOBRA2 59 (13) 
HsCOBRA3 60 (12) 
HsCOBRA4 59 (13) 
HsCOBRAS 295(59) 
HsCOBRA6 92(75) 
HsCOBRA7 97 (81) 
HsCOBRA8 63(44) 
HsCOBRA9 11(8) 
HsCOBRA10 14(17) 
HsCOBRA11 101(40) 
HsCOBRA12 35(23) 
HsCOBRA13 96(28) 
HsCOBRA 14 59(9) 
HsCOBRAI5 34(45) 
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100 
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EH 19 种 三 级 结构 模型 ，HsCOBRA7 氨基 酸 比例 最 高 (81%) , HsCOBRA13, 
HsCOBRA24 均 具 有 水 解 酶 结构 ，HsCOBRA1、HsCOBRA6、 
| 剂 结 构 , HSCOBRA2. HsCOBRA3., HsCOBRA4, HsCOBRA14 
有 碳水 化 合 物 结构 。HsCOBRA6、HsCOBRA7、HsCOBRA17 序列 同 源 性 为 99.9%; 


E 物 学 功能 ; HsCOBRAI. 
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图 3 蛇 足 石 杉 COBRA 家 族 三 级 结构 预测 图 
Fig.3 The tertiary structure prediction of COBRA family of Huperzia serrata 


2.6 te EA COBRA 基因 家 族 顺 式 作 用 元 件 分 析 


顺 式 作 


rH 
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4) , HPA 23 和 有 


牛 分 析 结 果 显 示 ，HsCOBRAs 共有 46 fh 
I EH 


具有 功能 注释 的 顺 式 作 用 元 件 ( 图 
Hc t+ BS JEM PY (MRE, Box 4, G-box, ACE, circadian, AE-box, 


GTI-motif, Spl、 3-AF1 binding site, CAG-motif. GA-motif. Gap-box. TCT-motif. chs-CMA2a, 


chs-CMA2c 、 
GTGGC-motif, TCCC-motif) 。 部 分 蛋 

CGCN4 motif) 及 响应 脱落 酸 (ABRE)、 生 长 素 (TGA-elem 
水 杨 酸 (TCA-element) 元 件 ， 此 外 还 含有 响应 
早 (MBS) 、 低 温 (LTR) 等 环境 的 顺 式 作 


为 局 动 子 和 增强 子 


Z 


GATA-motif ~ LAMP-element 、 


的 顺 式 作 


[X Jak 


元 人 


Box II, I-box., ATl-motif、 ATCT-motif、 
有 参与 分 生 组 织 表 达 (CAT-box) 和 胚乳 表达 


ent), JÆ A (P-box,. TATC-box ), 
大 氧 CARE) 、 缺 氧 (GC-motif)、 


Fo CAAT-box 元 件 的 数量 最 多 (307 个 )， 


有 的 顺 式 作 用 元 件 , 分 布 于 所 有 和 蛋白; 其 次 是 作为 核心 启动 顺 式 作 


用 元 件 的 TATA-box (206 个 ) 和 参与 干旱 诱导 MYB 结合 位 点 的 MBS (72 A) . HSCOBRAI8 
所 含 的 顺 式 作 用 元 件 最 多 (78 个 ) ， 其 次 为 HsSCOBRA14 (75 个 ) 和 HsCOBRA24 (70 个 )， 
HsCOBRA7 所 含 的 顺 式 作用 元 件 最 少 C10 个) 。 部 分 顺 式 作用 元 件 仅 在 一 条 蛋白 中 显示 ， 
如 TCCC-motif ( HsCOBRAS ) , 3-AF1 binding site ( H3COBRA12) , GTGGC-motif 
CHsCOBRA17) , ACA-motif, Box II -like sequence (HsCOBRAI8) , Box II, chs-CMA2c 
CHsCOBRA19) , ATCT-motif. ATl-motif. A-box (HsCOBRA20) , GA-motif. ACE 
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Partial cis-acting elements overlap. 


4 WEERT COBRA 家 族 顺 式 作用 元 件 预测 图 
Fig. 4 Cis-acting elements prediction of COBRA family of Huperzia serrata 


2.7 HEE AA COBRA 基因 家 族 保守 基 序 分 析 
HsCOBRAs 保守 基 序 结果 表明 ， 该 家 族 共 包含 6 种 保守 基 序 ， 命 名 为 motif I-motif 6 
(图 $， 图 6) 。 其 中 ， 仅 有 13 个 蛋白 拥有 motif，3 个 蛋白 (HsCOBRA2, HsCOBRA3, 
HsCOBRA4) 所 含 的 motif 最 多 , 均 为 6 种 ;7 个 蛋白 (HsSCOBRA1、HsCOBRA5、HsCOBRA7、 


HsCOBRA9, HsCOBRA10, HsCOBRA12. HsCOBRA24) 仅 有 1 种 motif. motif 6 分 布 最 
广 ， 其 次 为 motif 1. motif 3. motif 4. motif 3 含有 高 度 保守 的 CCVS 44, 7 个 蛋白 具有 
motif 3. HsCOBRAl. HsCOBRA6. HsCOBRA7 (motif 4) , HsCOBRA2. HsCOBRA3. 
HsCOBRAA (motif 1~motif 6) , HsCOBRA9, HsCOBRA10 (motif 6) 分 别 具 有 相同 种 类 的 
motif. 
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图 5 蛇 足 石 杉 COBRA 家 族 保守 基 序 预测 图 
Fig. 5 Conserved motif prediction of COBRA family of Huperzia serrata 


motif 4 


mofi 5 


图 6 WERI COBRA 家 族 保守 基 序 序列 分 析 


Fig.6 Conserved motif sequence analysis of COBRA family of Huperzia serrata 


2.8 ERT COBRA 基因 家 族 组 织 表达 模式 分 析 
根据 COBRA 基因 在 蛇 足 石 杉 叶 、 和 孢子 、 茎 、 芽 胞 中 的 转录 组 数据 分 析 其 表达 模式 ， 见 

图 7。HsCOBRA2 在 4 个 组 织 中 的 表达 量 均 最 高 ，HsCOBRA4 EH. fü. APN RIA 
次 于 HsCOBRA42。HsCOBRA47、HsCOBRA15、HsCOBRA20 在 4 个 组 织 中 几乎 不 表达 ， 此 外 
还 有 3 个 基因 (HsCOBRA6、HsCOBRA12、HsCOBRA23) 在 叶片 中 几乎 不 表达 ，1 个 基因 

(HsCOBRA21) 在 孢子 中 几乎 不 表达 ,4 个 基因 (HsCOBRA6、HsCOBRA12、HsCOBRA13、 
HsCOBRA23) 在 芭 中 几乎 不 表达 ，4 个 基因 (HsCOBRA4、HsCOBRA6、HsCOBRA21、 
HsCOBRA23) 在 芽 胞 中 几乎 不 表达 。COBRA 基因 家 族 在 蛇 足 石 杉 的 中 的 表达 量 为 叶 二 孢子 
> 茎 > 芽 胞 ; HsCOBRA2 和 HsCOBRA16 在 芽 胞 及 孢子 中 表达 量 较 高 ， 可 能 参与 昵 足 石 杉 的 
繁殖 。 


HcoBRt 
加 3409| 15.01| 0.00 | HsCOBRA4 180.00 
| 1.58 | 


IN: HsCOBRA3 150.00 


HsCOBRAM 


| 0.96 ES7771E Tx 
Lose [023 | 032 | 023 | HsCOBRA24 dia 
HsCOBRAY9 3000 


Foo oi oo | ooo We S 
Fooo | 002 [090 | 0.00 BERE 
Fooo | 000 [000 | 0.00 Wt 
Fooo | 0.00 0.00 | 0.00 We 
HSCOBRAVI 
1.86 HsCOBRAS 
0.98 HsCOBRA\0 
HsCOBRAI 
| 4.28 | 2.37 -cosrs 
HsCOBRA14 


HsCOBRAIG 
A c D 


B 


A. Ut; B. fü; C. 4; D. 芽 胞 。 
A. Leaves; B. Spores; C. Stems; D. Gemma. 
图 7 蛇 足 石 杉 COBRA 家 族 的 组 织 表达 模式 分 析 
Fig. 7 Expression profiles of COBRA family of Huperzia serrata 


2.9 Re EAI COBRA 基因 家 族 系统 发 育 分 析 
为 了 解 蛇 足 石 杉 与 蕨 类 植物 COBRA 基因 家 族 的 进化 关系 ， 构 建 蛇 足 石 杉 与 江南 卷 柏 
(Selaginella moellendorffii) 的 COBRA 家 族 成 员 的 NJ 系统 发 育 树 ， 见 图 8。 序列 比 对 发 现 
HsCOBRA2、HsCOBRA3、HsCOBRA4、HsCOBRA11、HsCOBRA13、HsCOBRA14 具有 
CCVS 结构 ， 分 别 位 于 蛋白 序列 的 第 238 位 和 第 241 位 、 第 268 位 和 第 271 位 、 第 238 位 和 
第 241 位 、 第 20 位 和 第 23 位 、 第 109 位 和 第 112 位 、 第 431 位 和 第 434 位 氨基 酸 残 基 中 ， 
蛋白 高 度 保 守 。 图 8 结果 表明 ，HsCOBRA1 与 HsCOBRA7, HsCOBRA2 与 HsCOBRA4, 


HsCOBRAI16 与 HsCOBRAI7 具有 同 源 性 ， 蛇 足 石 杉 与 江南 卷 柏 COBRA 基因 家 族 成 员 之 间 
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图 8 蛇 足 石 杉 与 江南 卷 柏 COBRA 家 族 的 系统 发 育 树 
Fig. 8 Phylogenetic tree of COBRA family of Huperzia serrata and Selaginella moellendorffii 


3 讨论 与 结论 


ACHE FFE HE LH 24 SE Ee COBRA 基因 家 族 成 员 ， 且 多 以 碱 性 蛋白 为 主 ， 二 级 结 
构 主 要 由 无 规则 卷曲 和 a 螺旋 及 延长 链 构 成 ， 大 部 分 成 员 具 有 交响 应 顺 式 作用 元 件 ， 与 大 白 
菜 ( 廉 瑞 婷 等 ，2022) 、 亚 麻 ( 齐 燕 妮 等 ，2019)、 高 梁 〈 元 志 成 等 ，2020) 、 番 茄 ( 曹 颖 等 ， 
2012)、 玉 米 ( 潘 倚天 ， 黄 敏 ，2022)、 茶 树 ( 艾 安 涛 等 ，2021) 等 COBRA 基因 家 族 相 似 ， 暗 示 
这 些 物种 中 可 能 具有 数量 相近 、 人 性 质 和 功能 相似 的 家 族 成 员 。 但 研究 也 发 现 HsCOBRAs XE 
细胞 定位 与 玉米 、 亚 麻 、 谷 子 ( 任 昂 彦 等 ，2021) 中 COBRA 定位 不 同 ， 提 示 不 同 物种 的 基因 
在 进化 过 程 中 可 能 产生 了 某 种 特定 功能 的 分 化 。HsCOBRAs 中 6 个 蛋白 (HsCOBRAIS, 
HsCOBRA19, HsCOBRA20, HsCOBRA21, HsCOBRA22, HsCOBRA23) 无 结构 域 及 保守 
基 序 , 且 在 蛇 足 石 杉 四 个 组 织 中 几乎 不 表达 , 该 6 个 蛋白 可 能 不 属于 蛇 足 石 杉 COBRA 家 族 


e 


成 员 。 
石 杉 碱 甲 在 治疗 阿尔 茨 海 默 病 及 益 智 中 表现 出 时 效 长 、 副 作用 少 的 优势 〈 陈 思 思 等 ， 
2021) ， 蛇 足 石 杉 作为 石 杉 碱 甲 的 主要 来 源 之 一 ， 其 原 植物 的 生长 繁殖 、 产 量 与 可 持续 发 展 


备 受 人 们 关注 。COBRA 基因 主要 功能 是 参与 植物 细胞 扩 增 ( 廉 瑞 婷 等 ， 2022) ， 通 过 调节 植 
物 的 生长 发 育 从 而 影响 其 生物 量 和 产量 。 研 究 中 发 现 HsSCOBRAs 具有 参与 分 生 组 织 表 达 及 
胚乳 表达 的 顺 式 作 用 元 件 ， 且 16 个 HsCOBRAs EI, WF. Z, FENER RAA EE 


FENZIA, CHAM RAP PRA, UE ea COBRA 基因 家 族 参与 到 其 叶 、 抱 
子 、 茎 及 芽 胞 的 发 育 过 程 ( 杨 兰 锋 等 ，2023) 。 前 期 研究 发 现 ， 蛇 足 石 杉 大 叶 、 老 叶 及 孢子 
圳 壁 中 均 具 有 丰富 的 纤维 素 ， 提 示 HsCOBRAs 基因 家 族 可 能 通过 调节 纤维 素 及 细胞 壁 生长 、 
机 械 强 度 等 参与 蛇 足 石 杉 的 生长 发 育 过 程 , 提高 其 生物 量 。 移 子 是 蛇 足 石 杉 主 要 繁殖 器 官 之 
—, HsCOBRAs 基因 家 族 可 能 影响 其 爷 子 宫 壁 的 生长 发 育 从 而 影响 着 孢子 的 成 熟 和 萌发 。 同 
时 ， 该 基因 家 族 成 员 具 有 响应 干旱 、 缺 氧 、 低 温 、 盐 碱 等 非 生 物 胁迫 的 顺 式 作用 元 件 。 其 中 
ACA 能 响应 冷 、 热 、 盐 胁迫 和 紫外 等 多 种 非 生物 胁迫 《〈 陈 逃 飞 等 ，2023 ) ; ABRE 元 件 能 
a a NON LUN 
LTR Spee Pc AEF BE UI Ek ABA, H20» 和 和 干旱 等 而 发 生 转 录 水 平 的 变化 ( 周 宾 寒 等 ,2023 ) 。 
ACA. ABRE, LTR 反 转 录 转 座 子 在 HsCOBRAs 基因 家 族 启动 子 区 域 大 量 存在 ， 说 明 蛇 足 石 
Mauer Ae awa pointe UE atte ie IR 

蛇 足 石 杉 COBRA 基因 家 族 具 有 与 其 他 物种 相似 的 生物 信息 学 特征 及 规律 ， 同 时 
能 参与 生长 发 育 、 繁 殖 及 响应 不 良 环境 胁迫 的 生物 学 功能 。 通 过 下 Fite LAG COBRA 基 
家 族 的 分 子 生物 信息 学 特征 及 组 织 表 达 规 律 ， 有 助 于 阐明 蛇 足 石 杉 COBRA 基因 的 分 子 调控 
和 促进 生长 发 育 机 制 ， 为 生物 量 的 增长 提供 新 途径 ,为 将 来 的 分 子 育 种 及 基因 功能 验证 等 和 
完 提 供 理 论 依据 。 
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